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Present Status of the Muon Anomalous Magnetic Moment 

Eduardo de RafaeP* 

a Centre de Physique Theorique, CNRS-Luminy, Case 907, F-13288 Marseille Cedex 9, France 

These pages, based on my talk at the Montpellier 14th International Conference in QCD, provide us with a 
short update of the Standard Model contributions to the muon anomalous magnetic moment. 
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1. Introduction 

We shall be concerned with the g-factor of the 
muon which relates its spin to its magnetic mo- 
ment 



2m u c 




(1) 



more precisely with the correction a M to the Dirac 
value — 2 which generates the so called anoma- 
lous magnetic moment. The present experimental 
world average determination, which is dominated 
by the latest BNL experiment (the E821 collabo- 
ration [1]), is 

a° xp = 116 592 080(63) x 10~ n (0.54 ppm) , (2) 

where the origin of the error is 0.46 ppm statisti- 
cal and 0.28 ppm systematic. This determination 
assumes CPT-invariance i.e., a„- = a^+. 

The question we shall discuss is: how well can 
the Standard Model digest this precise number?. 
As we shall see, the precision of a° xp is such that 
it is sensitive to the three couplings of the Gauge 
Theory which defines the Standard Model, as well 
as to its full particle content d. 

2. The QED Contributions (Leptons) 

This is by far the dominant contribution, which 
is generated by two types of Feynman diagrams: 
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ref. [2]. 



2.1. The Massless Class 

This class is generated by Feynman diagrams 
with virtual photons only as well as by diagrams 
with virtual photons and fermion loops of the 
same flavour as the external particle (the muon 
in our case) . Since a M is a dimensionless quantity, 
this class of diagrams gives rise to the same con- 
tribution to the muon, electron and tau anoma- 
lies. They correspond to the entries a^ 2 ") in Ta- 
ble 1, with n— 1, 2, 3, 4 indicating the number of 
loops involved. They are known analytically at 
one loop [3]; two loops |4|5| : and three loops [6]. 
This is the reason why there is no error in the 
corresponding numbers in the second column of 
Table 1. 

At the four-loop level, there are 891 Feynman 
diagrams of this type. Some of them are already 
known analytically, but in general one has to re- 
sort to numerical methods for a complete evalu- 
ation. This impressive calculation, which is sys- 
tematically pursued by Kinoshita and collabora- 
tors, requires many technical skills and is under 
constant updating; in particular thanks to the ad- 
vances in computing technology. The entry 
in Table 1 is the one corresponding to the most 
recent published value [7], with the error due to 
the present numerical uncertainties. 

Notice the alternating sign of the results from 
the contributions of one loop to four loops, a sim- 
ple feature which is not yet a priori understood. 
Also, the fact that the sizes of the (f)" coeffi- 
cients for n = 1,2, 3,4 remain rather small is an 
interesting feature, allowing one to expect that 
the order of magnitude of the five-loop contribu- 
tion, from a total of 12 672 Feynman diagrams, 
is likely to be of 0(a/-K) 5 ~ 7 x 10~ 14 . This 
is well beyond the accuracy required to compare 
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Table 1 

QED Contributions (Leptons) {a" 1 = 137.035 999 084 (51) [0.37 ppb] [8]} 



Contribution 


Result in Powers of - 


i 


Numerical Value in 10 11 Units 


a& 


0.500 000 000 (00) (f 


) 


116 140 973.29 (0.04) 


aW 


-0.328 478 965 (00) (f ) 


'2 

1 




a^ 4) (total) 


0.765 857 410 (27) (f ) 


2 

1 


413 217.62 (0.01) 


W 


1.18 124 146 (00) (f ) 


, 3 




4 6) (total) 


24.05 050 964 (43) (f ) 


3 

1 


30 141.90 (0.00) 


a& 


-1.9 144 (35) (f) 


4 

1 




(total) 


130.8 055 (80) (f ) 


4 

1 


381.33 (0.02) 


a\} ^ (total estimate) 


663 (20) (f ) 




4.48 (1.35) 


a (2 + 4 +6+8+ 10) (QED) 






116 584 718.09 (0.14)(0.04) 



with the present experimental result for a M , but 
it will be eventually needed for a more precise 
determination of the fine-structure constant a 
from the precision measurement of the electron 
anomaly |8J. 

2.2. The Massive Class 

This second class is generated by Feynman di- 
agrams with lepton loops of a different flavour to 
the one of the external muon line. Their con- 
tribution to is then a function of the lepton 
mass ratios involved. These contributions are 
generated by vacuum polarization subgraphs and 
by light-by-light scattering subgraphs involving 
electron and tau loops. Both the two-loop and 
three-loop contributions of this class are known 
analytically d. The full three-loop evaluation in- 
volving electron-loop subgraphs, by Laporta and 
Remiddi |9I10| . is a remarquable achievement. 
The numerical errors quoted in Table Q] for these 
contributions are due to the present experimental 
errors in the lepton masses |llj . 

At the four-loop level, only a few contributions 
are known analytically. Kinoshita and his col- 
laborators have, however, accomplished a full nu- 
merical evaluation of this class (see ref. [12J and 
references therein.). The corresponding error in 
Table[T]is the combined error in the lepton masses 
and the present error due to the numerical inte- 
gration. 

3 For a history of the successive improvements in the eval- 
uation of these contributions see e.g. ref. [2]. 



The number quoted for the full five-loop QED 
contribution in Table Q] is the present estimate 
quoted in ref. [I3j. It is likely to be improved in 
the near future. 

2.3. The Mellin Barnes Technique 

There has been a recent technical development 
in the evaluation of Feynman diagrams involving 
mass ratios, which has already been useful in the 
evaluation of some higher order contributions to 
cifj_ (see refs. |14|15j ) and which seems promis- 
ing for further calculations. In these papers it is 
shown how the Mellin-Barnes integral represen- 
tation of Feynman parametric integrals allows for 
an easy evaluation of as many terms as wanted in 
the asymptotic expansion of Feynman diagrams 
in terms of one and two mass ratios. 

The basic idea is to express the contribution 
to a M from a Feynman diagram, or a class of di- 
agrams, as an inverse Mellin transform with re- 
spect to the mass ratios involved in the diagrams. 
The remarkable property of this representation is 
the factorization in terms of moment integrals. It 
is in fact this factorization which is at the basis 
of the renormalization group properties discussed 
in ref. [16], and used since then by many authors 
(see e.g. ref. [17] and references therein). The al- 
gebraic factorization in the Mellin-Barnes repre- 
sentation, however, is more general. The stan- 
dard renormalization group constraints only ap- 
ply to the evaluation of asymptotic behaviours in 
terms of powers of logarithms and constant terms. 




Fig.l Diagrams with three e-loops and a r-loop. 

In the Mellin-Barnes framework, this appears as 
a property of the residue of the leading Mellin 
singularity. What is new here is that this ex- 
tends as well to the subleading terms, which are 
governed by the residues of the successive Mellin 
singularities (in the negative real axis, in the case 
of electron loops) ; or by two-dimensional residue 
forms |15|18| . in the case of the Mellin singular- 
ities associated to two mass ratios (the case of 
both electron and tau loops). 

As an example, we quote a few terms of the 
result obtained for the tenth-order contribution 
from the string of vacuum polarization subgraphs 
shown in Fig. 1: 
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In fact, the analytic calculation in ref. [H] which 
leads to this precise number, also includes terms 



up to O 



(2 \ *± 2 
— £ log J -f 



which are already 



smaller than the error generated by the lepton 
masses in the leading order terms. 

3. Hadronic Contributions 

The electromagnetic interactions of hadrons 
produce contributions to induced by the 
hadronic vacuum polarization and by the 
hadronic light-by-light scattering. 



Hadrons 



Fig. 2 Hadronic Vacuum Polarization 

3.1. Hadronic Vacuum Polarization 

All calculations of the lowest-order hadronic 
vacuum polarization contribution to the muon 
anomaly (see Fig. 2) are based on the spectral 
representation [19| 



hvp 



a r°° dt i , x r 1 

- / Imn(t) / dx- 

TT Jo t TT J a 



x 2 (l - x) 
■ + -V(l-x) 



(4) 



with the hadronic spectral function ilmll(t) re- 
lated to the one-photon e + e~ annihilation cross- 
section into hadrons (m e — > 0) as follows: 



A-K 2 a 1 



* (7) — > hadrons} 



t 



Imll(t) . 



(5) 



This contribution is dominated by the 7r + 7r~ 
channel; the region of the p-resonance in parti- 
cular |20|21j . The history of evaluations of ajf p 
is a long one which can be traced back, e.g. in 
ref. [2]. The most recent compilation of e + e~ an- 
nihilation data used in the evaluation of the dis- 
persive integral in Eq. (@]) made by Michel Davier 
and collaborators, which also includes the new 
precise measurements from the experiments SND 
and CMD-2 at Nobosibirsk as well as some exclu- 
sive channels from BaBar, gives the result Q: 



7 hvp _ 



(6 908±39 eX p±19 rad ±7 Q cD) x 10 



-11 



(6) 



Unfortunately, the discrepancy with the evalua- 
tion made using the r-spectral functions, cor- 
rected for isospin-breaking effects, still persists. 
Here, one has to wait for the forthcoming re- 
sults from the high precision measurements on 
the 7T7T mode at BaBar using the radiative return 
method. We shall then be able to check the con- 
sistency with the result in Eq. ([6]) and, hopefully, 
improve the accuracy. 



See e.g. ref [22] and references therein for details. 
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There is a similar spectral representation to 
the one in Eq. (@]) for the next-to-leading order 
hadronic vacuum polarization [23], with the ker- 
nel |24|25] in Eq. (j4]) , replaced by a two-loop ker- 
nel, which is also known analytically [26]. The 
most recent numerical evaluation, using the same 
data as for the lowest-order evaluation, gives 



a hv P (nio) = (_ 97 .9 ± o.9 exp ±0.3 r 



rad 



ixlO" 



(7) 



3.2. Hadronic Light— by— Light Scattering 

Unlike the hadronic vacuum polarization con- 
tribution, there is no direct experimental input 
for the hadronic light-by-light scattering contri- 
bution to shown in Fig. 3; therefore one has 
to rely on theoretical approaches. 

So far, the only rigorous theoretical result is the 
observation that, in the QCD large-iV c limit and 
to leading order in the chiral expansion, the dom- 
inant contribution comes from the Goldstone-like 
neutral pion exchange which produces a characte- 
ristic universal double logarithmic behavior with 
a coefficient which can be calculated exactly [27J: 



K ' Vtt/ 48tt 2 F2 L m. 



L + 0(ln?*)+0®\ (8) 



where F v denotes the pion coupling constant in 
the chiral limit {F„ ~ 90 MeV). Testing this limit 
was particularly useful in fixing the sign of the 
phenomenological calculations of the neutral pion 
exchange [28]. 

Although the coefficient of the In (m p /m w ) 
term in Eq. ([8]) is unambiguous, the coefficient 
of the ln(m p /m7r) term depends on low-energy 
constants which are difficult to extract from ex- 
periment |27|29] (they require a detailed know- 
ledge of the 7T° — > e + e" decay rate with inclusion 
of radiative corrections). Moreover, the constant 
term in Eq. ([5]) is not fixed by chiral symmetry re- 
quirements, which makes the predictive power of 
an effective chiral perturbation theory approach 
rather limited for our purposes. Therefore, one 
has to adopt a dynamical framework which takes 
into account explicitly the heavier meson degrees 
of freedom as well. 

The most recent calculations of a^ 11 in the 
literature [28 30 31 32] are all compatible with 
the QCD chiral constraints and large-iV c limit 




Fig. 3 Hadronic Light-by-Light Scattering 

discussed above. They all incorporate the 7r°- 
exchange contribution modulated by 7r°7*7* form 
factors, correctly normalized to the ir° — > 77 de- 
cay width. They differ, however, in the shape 
of the form factors, originating in different as- 
sumptions: vector meson dominance in a specific 
form of Hidden Gauge Symmetry in Ref . [30] ; in 
the form of the extended Nambu-Jona-Lasinio 
(ENJL) model in Ref. pffj; large-iV c models in 
Refs. |28|32] : and on whether or not they satisfy 
the particular operator product expansion con- 
straint discussed in Ref. [32]. 

In order to compare different results it is con- 
venient to separate the hadronic light-by-light 
contributions which are leading in the l/N c - 
expansion from the non- leading ones [33j . Among 
the leading contributions, the pseudoscalar me- 
son exchanges which incorporate the ir°, and to 
a lesser degree the r\ and rf exchanges, are the 
dominant ones. As discussed above, there are 
good QCD theoretical reasons for that. In spite 
of the different definitions of the pseudoscalar me- 
son exchanges and the associated choices of the 
form factors used in the various model calcula- 
tions, there is a reasonable agreement among the 
final results. The result quoted in a recent update 
discussed in ref. [34] gives: 



a hn (n ,r) ,r]') = (114 ±13) x 10" 



(9) 



Other contributions, which are also leading in the 
l/-/V c -expansion, due to axial-vector exchanges 
and scalar exchanges, give smaller contributions 
with updated errors, as discussed in ref. |34J: 



a hll (l+) = (15 ±10) x 10- 



and 



a 



hll (0+) = -(7 ±7) x 10" 



(10) 



(11) 



The subleading contributions in the 1/N C - 
expansion are dominated by the charged pion 
loop. However, because of the model dependence 
of the results one obtains when the pion loop is 
dressed with hadronic interactions it is suggested 
in ref. [34] to use the central value of the ENJL 
model evaluation in [31], but with a larger er- 
ror which also covers unaccounted loops of other 
mesons, : 

a hn (7T+n-) = -(19 ±19) x KT 11 . (12) 

From these considerations, adding the errors in 
quadrature, as well as the small charm contribu- 
tion: a hll (c) = 2.3 ± xlCT 11 , one gets 

a hu = (105 ±26) x 1CT 11 , (13) 

as a final estimate. 

4. Electroweak Contributions 

The leading contribution to a M from the Elec- 
troweak Lagrangian of the Standard Model, ori- 
ginates at the one-loop level. The relevant Feyn- 
man diagrams (in the unitary gauge) are shown 
in Fig. 4. The analytic evaluation of the overall 
contribution gives the result (see e.g. ref. |35j): 




= 194.8 x 1CT 11 , (14) 



where the weak mixing angle is defined by 
sm 2 6 w = 1 - M^/Mf ~ 0.223, and G F ~ 
1.166 x 10 -5 is the Fermi constant. Notice 
that the contribution from the Higgs boson, 
shown in parametric form in the second line, 

is of O f%?4-Si- In ^¥), rather small for the 

V y/2 4tt 2 Mfj m 2 ) ' 

present lower bound on Mh, but known analyti- 
cally. 

The a priori possibility that the two-loop 
electroweak corrections may bring in enhance- 
ment factors due to large logarithms, like 
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Fig. 4 Weak interactions at the one-loop level 



ln(Mf/m^) ~ 13.5, has motivated a thorough 
theoretical effort for their evaluation, which has 
been quite a remarkable achievement. 

It is convenient to separate the two-loop elec- 
troweak contributions into two sets: those con- 
taining closed fermion loops and the bosonic 
corrections, which we denote by (bos). 
The latter have been evaluated using asymp- 
totic techniques in a systematic expansion in 

2 M 2 M 2 

powers of sin 9w, where log — f- terms, log -j-M- 

M 2 M 2 M 2 

terms, -r-#- log -jit- terms, -y# terms, and con- 

stant terms are kept. Using sin 2 9w = 0.223 and 
50 GeV < M H < 700 GeV results in [36137)38] : 

a W) (bos) = ^|j| x - ( _ 82 . 2± 5.9) 
V2 071"* tt 

= (-22.2 ± 1.6) X 1CT 11 . (15) 

The discussion of the fermionic corrections is 
more delicate. Because of the U(l) anomaly can- 
cellation between lepton loops and quark loops 
in the electroweak theory, one cannot separate 
hadronic from leptonic effects any longer in di- 
agrams like the ones shown in Fig. 5, where a 
VVA-triangle with two vector currents and an 
axial-vector current appears. It is therefore ap- 
propriate to separate the fermionic corrections 
into two subclasses. One is the class in Fig. 5, 
which we denote by a^ W ^ (I, q) . The other class 
is defined by the rest of the diagrams, where 
quark loops and lepton loops can be treated sepa- 
rately, which we call a^ W ^ (ferm-rest). This latter 
contribution has been estimated to a very good 
approximation in ref. [36] with the result, 

£t„ '(ferm-rest) = — = — hr — X (-21 ± 4) , (16) 

where the error here is the one induced by dia- 
grams with Higgs propagators with an allowed 
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Table 2 Standard Model Contributions 



Fig. 5 Two-loop electroweak diagrams generated by 
the 77Z -Triangle. There are similar diagrams cor- 
responding to the /i, c, s and r, t, b generations. 

Higgs mass in the range 114 GeV < Mh < 
250 GeV. 

Concerning the contributions to af W ^ 2 \l, q), it 
is convenient to treat the three generations sep- 
arately. The contribution from the third genera- 
tion can be calculated in a straightforward way 
using effective field theory techniques [39], be- 
cause all the fermion masses in the triangle loop 
are large with respect to the muon mass, with the 
result |39|36| : 



I r (2) (-,M) = %^-x(-3o. 6) . 

V2 ott z n 



(17) 



However, as first emphasized in ref. [39], an ap- 
propriate QCD calculation when the quark in the 
loop of Fig. 5 is a light quark should take into ac- 
count the dominant effects of spontaneous chiral- 
symmetry breaking. Since this involves the u , d 
and s quarks, it is convenient to lump together 
the contributions from the first and second gen- 
erations. An evaluation of these contributions, 
which incorporates the QCD long-distance chi- 
ral realization |39|40] as well as perturbative [41] 
and non-perturbative 10117] short-distance con- 
straints, gives the result 



Ewm, , \ 

a n ( e > Mi u i "i s > c ) 



Op m£ a 
^2 8tt 2 tt 



x(-24.6±1.8). (18) 



From the theoretical point of view, this calcula- 
tion has revealed surprising properties concern- 
ing the non-anomalous component of the VVA- 
triangle [42], resulting in a new set of non- 
renormalization theorems in perturbation the- 
ory [42143] . 

Putting together the partial two-loop results 
discussed above, one finally obtains for the overall 
electroweak contribution the value 



Contribution 


Result in 10 11 units 


QED (leptons) 


11 6584 718.09 ± 0.14 ± 0.04 Q 


HVP(lo) 


6 908 ± 39 exp ± 19„ d ± 7 pQCD 


HVP(ho) 


-97.9±0.9 exp ±0.3 rad 


HLxL 


105 ± 26 


EW 


152 ± 2 ± 1 


Total SM 


116 591 785 ± 51 



EW(l) , Gf m u 

152(2)(1) x 10 -11 . 



-158.4(7.1)(1.8)] 



(19) 



where the first error is essentially due to the Higgs 
mass uncertainty, while the second comes from 
hadronic uncertainties in the VVA-loop evalua- 
tion. The overall result shows indeed that the 
two-loop correction represents a sizeable reduc- 
tion of the one-loop result by an amount of 22%. 
An evaluation of the electroweak three-loop lead- 
ing terms of O 



In 



M z 



malization group arguments [4414 



, using renor- 

], shows that 
higher order effects are negligible [0(1O -12 )] for 
the accuracy needed at present. 

5. Summary 

Table 2 collects the various Standard Model 
contributions to a M which we have discussed. No- 
tice that the largest error at present is the one 
from the lowest order hadronic vacuum polariza- 
tion contribution. Adding experimental and the- 
oretical errors in quadrature gives a total 



,SM 



(116 591 785 ±51) x 10 



-li 



(20) 



with an overall error slightly smaller than the one 
in the experimental determination in Eq. Q • The 
comparison between these two numbers, shows an 
intriguing 3.6 a discrepancy. 
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